Disorder-Induced Vibrational Localization 
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The vibrational equivalent of the Anderson tight-binding Hamiltonian has been studied, with 
particular focus on the properties of the eigenstates at the transition from extended to localized 
states. The critical energy has been found approximately for several degrees of force-constant 
disorder using system-size scaling of the multifractal spectra of the eigenmodes, and the spectrum at 
which there is no system-size dependence has been obtained. This is shown to be in good agreement 
with the critical spectrum for the electronic problem, which has been derived both numerically and 
by analytic means. Universality of the critical states is therefore suggested also to hold for the 
vibrational problem. 

PACS numbers: 63.50.+X, 63.20.Pw 
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The Anderson electron localization problem Q| is one 
that has attracted much attention over the last 40 years. 
The fact that the problem can be stated so simply, and 
yet have startlingly complex consequences, has made it 
a challenging topic to work on (2). Indeed, only recently 
has it been possible to verify numerically many of the 
theoretical results on powerful supercomputers || . How- 
ever, the closely related vibrational problem has not been 
explored to the same degree, despite being similar enough 
to use the same techniques yet different enough to pro- 
duce new and interesting results. 

The phenomenon of localization is a second-order 
phase transition between eigenstates that are spatially 
localized and those that are delocalized, or extended |Q. 
In the thermodynamic limit, extended eigenmodes would 
cover the whole of space whereas localized eigenstates 
are those which only involve a local subset of the system 
within a typical localization length. In the crystalline 
case for both the electronic and vibrational problems, 
the eigenstates are simple Bloch states due to transla- 
tional invarience, and are therefore extended. For the 
electronic problem, disorder is generally introduced ei- 
ther in the on-site energy terms (diagonal disorder) or 
the interaction terms (off-diagonal disorder) In a 3d 
lattice with weak diagonal disorder, there are two critical 
energies, at the top and bottom of the band, at which 
the Localization-Delocalization (LD) transition occurs. 
As the degree of disorder is increased, these two critical 
energies approach, and finally meet. At this point, all 
the eigenmodes are localized and the system becomes an 
electrical insulator. This transition is termed the Metal- 
Insulator Transition (MIT) H. Off-diagonal disorder 
produces fundamentally different behaviour: at no level 
of disorder are all the eigenmodes localized, and hence 
there is no MIT §. 

Our approach to the problem of vibrational localiza- 
tion has been numerical, applying high-pcrfomance com- 
puters to the task of obtaining the eigenmodes. The 
Anderson electron Hamiltonian can be expressed in a 
site basis, giving a sparse matrix representation of the 



problem, for which the eigenvectors can then be found 
by using standard Lanczos methods. Modern computers 
can solve such eigenproblems for many millions of atomic 
sites. A bigger problem is how to recognise quantitatively 
the difference between localized and extended states. 

There are several methods for distinguishing extended 
from localized states, e.g. by looking at the properties 
of the Hamiltonian, such as the transfer matrix method 
H observing differences in the level-spacing statis- 
tics the Thoulcss criterion |{|, or by looking at the 
eigenstates themselves. The latter is not trivial though, 
since as the critical energy is approached from the local- 
ized regime, the localization length diverges. Thus, for 
a finite system size, the eigenmodes quickly become ex- 
tended over a larger range than the system size and it 
becomes difficult to assess whether a state is truly local- 
ized or extended. These states are known as prelocalized 
states [|l0|, and to characterize these as localized or ex- 
tended, we can use multifractal analysis (MFA) 

It has been suggested that the eigenmode at exactly 
the LD critical energy will show multifractal character- 
istics fl2|| . The standard way of characterising the mul- 
tifractality is the singularity spectrum, which has been 
shown for electrons |13| to be universal for an isotropic 
system (see Ref.jlJJ for treatment of an anisotropic sys- 
tem) and independent of the probability distribution of 
the disorder. The analytic predictions for the singularity 
spectrum pa] , based on the d = 2 + e expansion of the 
non-linear a model, are in good agreement with numerics 
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The aim of this paper is two-fold. Firstly, we use 
MFA in order to identify the threshold energy of the LD 
transition for different degrees of force-constant disorder 
and thus obtain the "phase diagram" in the frequency- 
disorder plane for vibrational excitations in disordered 
models. Secondly, we demonstrate the universal features 
of the multifractal critical states at the LD transitition 
for the vibrational problem. 

We can use the idea of critical multifractality to de- 
termine whether the states are extended or localized by 
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looking at how the singularity spectrum, characterising 
the eigenmode, changes with simulation-system size. For 
a true multifractal state, assuming that finite-size effects 
are small, the singularity spectrum will not depend on 
the simulation box size, whereas the spectra for states 
on either side of the LD transition will vary. Hence, by 
calculating the singularity spectrum for different system 
sizes, we can locate the critical energy [ fl3| . 

The harmonic vibrational problem that is addressed 
in this Letter can be formulated in a very similar way to 
the Anderson electron problem For vibrations, the 
equivalent to the electronic Hamiltonian is the symmetric 
dynamical operator: 
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(\i,a)-\j,a))((i,/3\-(j,l3\) , 



(1) 



with \i, a) being the site basis describing the displace- 
ment of atom i (i = 1 . . . N) along the Cartesian direction 
a (a — 1 . . . d, with d the dimensionality) . The matrix 
elements D 
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terms of force constants Kij, and unit vectors, Tij, 



Kij/2) (fij) a {fij)p are defined in 

con- 
necting the atoms i and j (for simplicity, all masses are 
taken to be equal, mj = 1). The dynamical matrix con- 
sists of d x d blocks with strong lattice symmetry-dictated 
correlations inside the blocks. Additionally, all the ele- 
ments of the on-diagonal blocks are the sums (with op- 
posite sign) of the similar elements of off-diagonal blocks, 
reflecting the sum-rule correlations in the dynamical ma- 
trix. Therefore, in the case of nearest-neighbor interac- 
tions considered below, the number of correlations be- 
tween the elements in the dynamical matrix is compara- 
ble with the number of independent random variables. 

There are three main differences between the Ander- 
son electron and the vibrational problems. Firstly, when 
there are no negative values of Kij, the system is me- 
chanically stable, so there are no negative eigenvalues, 
unlike the electron case. Secondly, the basic Anderson 
formulation gives a symmetric band structure. The vi- 
brational case is asymmetric for the third reason: there 
are d zero-frequency modes that cannot be localized since 
they correspond to bulk translational displacements of 
the system (Goldstone modes). Since the lower bound of 
the spectrum is therefore constrained to be extended in 
character, we expect in a single-band model that there 
will be only one LD transition near the high-frequency 
band edge. 

There are two major classes of model which can be 
used in studying localization: structures based on an un- 
derlying crystalline lattice with introduced disorder, and 
structures which have been created in an effort to recre- 
ate the distribution of atomic positions and bond angles 
found in real amorphous materials. For our study, we 
have chosen to analyse lattice models from the first class, 
with an underlying f.c.c. geometry and with the Kij in 



Eq. (|l|) taken from a probability distribution p(n). This 
is one of the simplest models and can be easily compared 
with the established results for the electron-localization 
problem for similar models. The distribution p{n) has 
been chosen to be a uniform (box) distribution, centered 
at /to = 1 with a full width A < 2kq in order to give both 
a simple random distribution and one where there are no 
negative force constants. Our models are face-centered 
cubic and range in size from L = 16 with 4096 atoms up 
to L = 48 with N = 110592 atoms. 

A multifractal is a generalization of a standard geo- 
metric fractal for the case when a single fractal dimen- 
sion cannot characterize the system [p.7| . For each point 
in our measure, we can define a value a(r) that describes 
the scaling of the measure with L around that point. 
If we now take the set of all points with a specific a, 
that itself is a fractal, with dimensionality f(a). The 
curve f(a) is known as the multifractal spectrum, or 
singularity spectrum, and can be used to characterize 
eigenmodes as localized or extended, as shown below. To 
calculate the scaling exponents, we define the measure 
Pi(Lb) = J2jeBo^i(L b ) l u j| 2 as tne sum of tlie squared 
displacements |uj| 2 of all the atoms j within the ith box 
of size Lb < L for a particular eigenmode, and examine 
how this measure scales with Lb, or equivalently, with 
A = Lb/L. We split our system up into N(X) boxes 
which completely and exactly cover the system, so that 
N(X) = \~ d . The standard normalization of the eigen- 
modes leads to a scaling law for the measure of the form 
(P(Lb)) L oc X d , averaging over all boxes. 

The assumption underlying multifractal analysis is 
that, for a finite interval of A, the qth moments of the 
P{L b ) also scale with power laws: (P q {L b )) L oc A d+r ^ 
where r(q) is independent of A. The range of A in our 
case has a lower bound at the interatomic spacing, since 
we are dealing with a discrete rather than a continuous 
system. The upper bound L/2 is dictated by finite-size 
effects. In the thermodynamic limit, as L — > oo (A — > 0), 
the states which satisfy the multifractal condition are 
only found exactly at the critical energy, and thus the 
exponents are defined uniquely as 



t(q) = lim 



In A 



(2) 



In practice, r(g) is found by performing a linear regres- 
sion of the calculated exponents with In A. From this we 
can obtain the singularity spectrum, f(a), where a is de- 
fined as a(q) = dr(q)/dq and f(a) is obtained from the 
Legendre transformation of r(g), f(a(q)) = a(q)q — T(q). 
Calculation of the singularity spectrum using the Legen- 
dre transformation suffers from numerical errors, so it is 
more convenient to calculate f(a) as a function of P(Lb) 
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explicitly |/i 



a(q) = 



— y 



P?(L b ) 



In A /-^ Z(q,L b ) 

boxes 



In A ^ Z(q,L b ) 

boxes 



lnPi(L 6 ), 



In P?(L 6 ), 



(3) 



(4) 



where Z(g,L 6 ) = £ bOTOS i?(L 6 ). 

Since we cannot take the limit A — > in Eqs.d|)-(||), 
the values of f(q) and a(q) are calculated by performing 
a linear regression of the respective sums with respect to 
In A. The linearity of these graphs is a good check of the 
multifractal nature of the measure. 

Care has to be taken over the box sizes used in the 
analysis. For example, taking the box to include just one 
atomic site proved to skew the regression, as did taking 
the box size be that of the entire system. The reason for 
the former is that the multifractality must break down 
at some point, certainly for box sizes on the order of 
the atomic spacing. Finite-size effects account for the 
discrepancy for the largest box size. 

The singularity spectra of the eigenmodes around the 
critical energy fluctuate strongly, and so it becomes nec- 
essary to take an average. Ideally, we would like to av- 
erage over different realizations of disorder, but in prac- 
tice this is only realistic for the smaller size models. For 
larger models, we take the computationally cheaper op- 
tion of averaging consecutive eigenmodes, which can be 
obtained simply in the Lanczos algorithm. In order to 
reduce errors, we have used the gliding-box method, av- 
eraging over all possible origins when dividing the system 
into boxes |Tij| . 

Once we have the spectra, we can find the frequency 
at which there is no change with system size to locate 
the mobility edge. Empirically, it was noted Jl4| that, 
for the Anderson plot of a(q) against (logL) 1 

gave a good linear fit with a different sign of gradient 
g(oj 2 ) = da(q;uj 2 )/d(logL)^ 1 on either side of the critical 
frequency. The same holds true for vibrational models, 
as clearly demonstrated in Fig. [|. We have therefore 
performed a linear regression on these curves, and the 
gradients of these lines have been plotted at different en- 
ergies to find the point where the singularity spectrum 
is size independent, at which g(uj 2 ) crosses the abscissa. 
We can get additional information by looking at different 
values of q. In practice, since a(q) is strongly correlated 
for similar q, we have looked at the representative values 
q = and 1, for which the g(oj 2 ) have opposite signs (see 
Fig. g. 

Initially, the analysis was undertaken throughout the 
acoustic band. We did not expect to find localization at 
the lower (zero-frequency) band edge Q, and indeed it 
was found that there was only one LD phase transition, 
located in the far high-energy band tail. The band edge 
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FIG. 1: Estimation of the localization edge for A = 1.5. Each 
line is at a different frequency, from uj 2 = 9.3 at the bottom 
to 9.5 at the top in steps of 0.02. The critical frequency is 
that at which this line has zero gradient. Notice the bold 
line shown, with approximately zero gradient, corresponds to 



9.4 and is at q(0) = 4.0 




FIG. 2: Plot of g = da(q)/d(logL)~ 1 for q = and 1. The 
squared critical frequency u> 2 is given by the zero-crossing 
point of the graph. In this case, u> 2 is between 9.4 and 9.44. 



calculated within the coherent potential approximation 
(CPA) was found to be quite close to the true localization 
threshold, as can be seen in Fig|], and therefore it can be 
used as a rough estimate for the frequency of the actual 
LD transition . 

Having found the mobility edge for several values of 
force-constant disorder A, we can plot these to produce a 
'phase diagram' of the eigenmodes. This is shown along- 
side the VDOS and the CPA band edge in Fig.g. As 
the localization edge is in the band tail and we are lim- 
ited to finite-size systems, few states are localized. With 
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- - CPA band edge 

— Mobility edge 

— - Crystalline VDOS 
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FIG. 3: Phase diagram showing the boundary between ex- 
tended and localized states. The VDOS for the crystal and 
the lattice with A = 2.0 have also been plotted to show the 
location of the mobility edge, ujI, within the band tail. The 
band edge calculated with CPA is also shown for reference. 




FIG. 4: Critical spectra for the force-constant disordered 
models. The parabolic approximation (PA) to Wegner's result 
[H is shown for comparison. 



increasing A, the mobility edge decreases in frequency 
with respect to the CPA band edge. However, since the 
band is broadening, the result is that the critical fre- 
quency actually increases with disorder, and thus there 
is no vibrational analogue to the electronic MIT. Simi- 
lar behaviour of the mobility edge with disorder can be 
seen in the phase diagram of the Anderson model with 
off-diagonal disorder H . 

For each degree of disorder, we obtain a new critical 
spectrum which is constant for each size. These critical 
spectra have been plotted in FigJ^, showing that for pos- 
itive values of q, i.e. the left hand side of the graph, all 



the spectra fit onto a master curve. The parabolic ap- 
proximation (PA) to Wegner's analytic result [Gl| is one 
which goes through the critical points f(a = 4) = 3 and 
f(a = 2) — 2, where the latter corresponds to the infor- 
mation dimension of the eigenmode. This PA has also 
been plotted on the graph for comparison. Note that the 
Wegner result is for the electronic Anderson model, yet 
it still fits well with the vibrational data, indicating a 
universality across the two different systems. The large 
error bars at high a are in the region where q is negative, 
where / and a are strongly dependent on the smallest 
values of the measure and where the errors in the eigen- 
modes themselves are largest. 

To conclude, we have investigated the localization phe- 
nomenon for vibrational excitations in disordered struc- 
tures, using an f.c.c. lattice model with force-constant 
disorder for analysis. Using MFA, we have confirmed 
the existence of only one LD transition in the acoustic 
band, and found the energy at which it occurs for differ- 
ent degrees of disorder. The eigenmodes at the threshold 
have been shown to be multifractal states exhibiting a 
quantitatively similar distribution function to that of the 
critical states in the electron Anderson model. 

We are grateful to R. Romer for supplying us with MFA 
code pofl , and to M.Schreiber for instructive communica- 
tions. 
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